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ABSTRACT 

We derive the Fundamental Plane (FP) relation for a sample of 1430 early-type galaxies in the 
optical (r band) and the near- infrared (K band), by combining SDSS and UKIDSS data. With such 
a large, homogeneous dataset, we are able to assess the dependence of the FP on the waveband. Our 
analysis indicates that the FP of luminous early- type galaxies is essentially waveband independent, 
with its coefficients increasing at most by 8% from the optical to the NIR. This finding fits well into 
a consistent picture where the tilt of the FP is not driven by stellar populations, but results from 
other effects, such as non-honiology. In this framework, the optical and NIR FPs require more massive 
galaxies to be slightly more metal rich than less massive ones, and to have highly synchronized ages, 
with an age variation per decade in mass smaller than a few percent. 
Subject headings: Galaxies: fundamental parameters - Galaxies: evolution 



1. INTRODUCTION 

Early-type galaxies (ETGs) populate a two- 
dimensional surface in the space of parameters 
that reflect size (effective radius), density (mean 
surface brightne ss), and kinetic temperature (veloc- 
ity dispersion) (iDjorgovskil Il992l). A key feature 



of the Fundam ental Plane (FP; jPressler et all Il987l : 
iDiorgovski fc D avis 19871) is its deviation ('tilt') from the 
virial theorem, which may be interpreted as a variation 
of the M/L ratio along the sequence of ETGs and/or 
the breaking of homology assumption, i.e. the fact that, 
for all galaxies, the observed parameters have the same 
power-law dependence on the corresponding physical 
quantities (namely, the central velocity dispersion on 
kinetic energy, the effective radius on 'gravitational' 
radius, and the effective surf ace brightness on the 
overall light profile, see e.g. IDiorgovski fc Santiagol 
[1993). Despite all the observational efforts, the origin 
of the tilt is still under debate. The change of the 
M/L ratio can be exp lained by a change in ei ther the 
stellar population (e.g. iPrugniel fc Simieiilll996D or dar k 
matter content with galaxy mass ( Ciotti et al.l Il996f ) . 
Both structural and dynamical non-homology have also 
been invoked as physical explaiiations of the observed 
tilt (s ee e.g. iHiorth fc MadsenI Il995l ICapelato et al " 
19951 iGraham fc CollessI Il99'7r an d iBusarello et"aL 



1997t l. Recently, iTrujillo et all (|2004 from now TBB04) 
showed that the tilt is mostly driven by dynamical and 
structural non-homology, while ste llar populations ac - 
count for only a small fraction of it. iBolton et al.l ()2007l ) 
argued that the tilt is more likely because of a variation 
of the dark matter content with mass, still favoring 
a picture where stellar populations play a minor role. 
Since the contribution of different stellar populations 
to galaxy luminosity is expected to be wavelength- 
dependent, while other effects (e.g. non-homology) are 
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not, the dependence of the FP on wavelength directly 
informs on how properties of the stellar populations 
change with mass, which is a crucial point to understand 
galaxy formation and evolution. 
Previous studies of the wa velength de p endence 

led to contradicto ry results. iPahre et al.l ()1998bl ) 

and lScodeggio et all ([Till) found the tilt to significantly 
decrease from optical to NIR wavelengths, interpreting 
this re sult as an increase of age and metallicity with 
mass. iMobasher eFan (|1999f ) and IZibetti et all (|2002l ) 
found only a small decrease of the tilt with wavelength, 
with the FP still being significan tly tilted in the NIR. 



However, iBernardi et al.l ( 2003b[ ). deriving the FP for 



ETGs observed in the Sloan Digital Sky Survey (SDSS), 
found evidence for the FP to be wavelength indepen- 
dent from the g to the z bands. Several different ef- 
fects can produce this puzzling picture. The FP in dif- 
ferent wavebands has often been derived for small sam- 
ples, with inhomogeneous measurements of galaxy pa- 
rameters, different selection criteria (e.g. galaxy samples 
spanning different ranges in magnitude and/or velocity 
dis persion), and with diff erent fitting methods. The FP 
bv IBernardi et all (|2003br ) avoided all these problems by 
analyzing the same sample of galaxies at different wave- 
bands, but it was limited to the short wavelength baseline 
provided by the SDSS. In the present work, for the first 
time, we derive the FP by using the same, large, homoge- 
neous sample of ETGs over the wide wavelength baseline 
provided by the r- and K-band data of the SDSS and the 
UKIRT Infrared Deep Sky Survey (UKIDSS). 

The layout of the paper is the following. In Sec. 2, we 
describe the selection of the sample, while Sec. 3 details 
how we obtain the r- and K-band structural parameters, 
and the central velocity dispersions. Sec. 4 deals with the 
comparison of the r- and K-band FPs. In Sec. 5, we show 
how the optical and NIR FPs constrain the variation of 
stellar population parameters along the galaxy sequence. 
The discussion follows in Sec. 6. Throughout the paper, 
we adopt the cosmology Hq — 75 km s^^ Mpc"-'^ , fim = 0.3, 
and riA = 0.7. 



2. SAMPLE SELECTION 
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We select a sample of ETGs, with available K-band 
photometry from the second data release of UKIDSS, and 
r-band photometry and central velocity dispersions from 
the fifth data release (DR5) of SDSS. First, a complete 
volume-limited catalog of galaxies is defined, consisting 
of all the 105036 objects in DR5 with an r-band abso- 
lute magnitude Mr< — 20, and a spectroscopic redshift 
in the range of 0.05 to 0.095. Absolute magnitudes are 
obtained from the SDSS r-band Petrosian magnitudes, 
k-correct ed to redshift 0.1 by using the fccorrectf 4_1_4 
software (jBlanton et al.ll2003[ ). The lower redshift limit 
is ch osen to minimize the aperture bias (jGomez et al.l 
|2003| ). while the upper r edshift limit guarantee s a high 
level of completeness (see ISorrentino et al.ll2006l) . ETGs 
are defined according to the SDSS spectroscopic parame- 
ter eclass, that classifies the spectral type on the basis of 
the principal component analysis technique, and the pho- 
tometric parameter fracDeVr, which measures the frac- 
tion of gala xy light that is fitted by a de Vaucouleurs law. 
Following iBernardi et~aI1 ()2003af ). we define as ETGs 
those objects with eclass <0 and fracDeVr >0.8, result- 
ing in a list of 47061 galaxies. Out of them, we retain 
only those 33628 galaxies with available central veloc- 
ity dispersion, ctq, between 70 and 420 kms~^. These 
cuts are required to obtain reliable a^'s from the SDSS- 
DR5 database. All the selected galaxies have spectra 
with median per-pixel S/N larger than 10, which is a 
further requirement to retrieve reliable SDSS velocity 
dispersions ^. The SDSS catalog is cross-matched with 
UKIDSS, resulting in 1570 galaxies. We notice that all 
the SDSS galaxies with reliable velocity dispersion, that 
are covered by the UKIDSS survey, are then included in 
this sample. In other words, the matching with UKIDSS 
does not change the magnitude limit of the present sam- 
ple. We select only those galaxies observed in K-band 
images with good seeing {FWHM < 1"). This selection 
reduces the sample size by only 10%, and excludes cases 
where structural parameters might be affected by large 
uncertainties. Since all the r-band images have FWHM 
smaller than 1.4" and the galaxies in our sample have 
on average effective radii larger in the optical than in the 
K-band (see below) , we do not apply any seeing selection 
to the SDSS photometry. The above procedure leads to 
a final sample of 1430 galaxies. 

3. FP PARAMETERS 

The photometric parameters entering the FP, namely, 
the effective radius, r^ and the mean surface bright- 
ness withi n that radius, < At > e; are derived using 
2DPH0T (jLa Barbera et al.ll2008f ). The r- and K-band 
images are processed by adopting the same 2DPH0T 
options, allowing homogeneous structural parameters to 
be derived between both bands. For each galaxy, a lo- 
cal PSF model is computed by fitting simultaneously the 
four closest stars to that galaxy on the image. Structural 
parameters, i.e. the eff'ective parameters and the Sersic 
index n (shape parameter), are then derived by modeling 
the galaxy images with two-dimensional seeing-convolved 
Sersic models. Effective radii are converted to physical 
units by using the angular diameter distance correspond- 
ing to the DR5 spectroscopic redshift, z, of each galaxy. 

^ See the list of requirements in the Algo- 
rithms section of the SDSS-DR5 website, at 

|http://www.sdss.org/dr5/algorithms/veldisp.html| 



Mean surface bri ghtnesses are de-reddened using the ex - 
tinction maps of iSchlegel. Finkbeiner and David (|1998f) . 
are corrected for cosmological dimming, by subtracting 
the term 7.51og(l-t-z), and are k-co rrected to redshift 0. 1 
with the fccorrectf 4_1_4 software (jBlanton et al.l 12003^. 
In order to estimate the accuracy on structural parame- 
ters, we use 160, out of the 1430, galaxies having repeated 
observations in different UKIDSS frames. We find that 
the averaged differences between repeated measurements 
of structural parameters are fully consistent with zero, 
amounting to -O.OliO.Ol, -0.037±0.04, and 0.005±0.01 
for logTc, </i>c and logn, respectively. The rms values 
of these differences amount to 32% in re, 0.6 mag/arcsec^ 
in </i>c, and 25% in n. Notice that the scatter in log re 
is fully consistent with the typical accuracy of the mea- 
surement of the half-light radii ( Kelson et al. 2000 ) . The 
quantity log Tg — 0.3 < fi >c, that is the relevant photo- 
metric parameter entering the FP, has an uncertainty 
of only 7%, as expected due to the correlation of mea- 
surement errors of the effective parameters. The com- 
parison of r- and K-band structural parameters for the 
present sample is fully con sistent with what is fou nd in 
our previous studies (e.g. iLa Barbera et al.l l2004t ). In 
particular, the mean ratio between r- and K-band radii 
is — O.llibO.Oldex, i.e. on average effective radii decrease 
by ^ 25% from the optical to the N IR. This valu e is con- 
sistent with that of -- 20% found bv lPahre et~an (lT998a), 
and is in agreement with the existence of negative color 
gradients in early-type galaxies. Sersic indices are fully 
consistent between optical and NIR. The average ratio of 
r- to K-band n values amounts to — 0.007ib0.009dex. Cen- 
tral velocity di spersions are retrieved from DR5 and are 
corrected as in IBernardi et al.] (| 2003br ) to a re lative aper- 
tureofr^/8, follow ing iJorgensen et al.l (|1995f l. As shown 
bylB erna rd1 (|2007D . for tJo < 150 km s"^ the DR5 velocity 
dispersions are slightly overestimated. This small bias in- 
creases up to 12% at (To ^ 100 km s~^. We remove this 
eff ect by applying the correction curve shown in fig. 4 
of iBernardil (|2007t ) (see the grey hue in the upper-left 
panel)^. The r- and K-band effective parameters, as well 
as the corrected DR5 velocity dispersions are given in 
Tab. 1 (fully available in electronic form) for all the 1430 
galaxies analyzed in the present study. Columns in the 
table provide the following quantities. Columns 1 and 2 
report right ascension (RA) and declination (DEC) in 
units of degrees. Columns 3 and 4 provide the logarithm 
of the effective radius (in unis of kpc) and the effective 
mean surface brightness (in unis of mag / arcsec^) in the 
r-band. Columns 5 and 6 list the same quantities as 
columns 3 and 4 but for the K-band. Column 7 reports 
the corrected DR5 velocity dispersions. 



4. THE SDSS AND UKIDSS FUNDAMENTAL PLANES 

We adopt the usual representation of the FP, log r-g — 
a log (To + b <^>c +c , where a and b are the "slopes", 
and c is the offset. These coefficients were derived by 
minimizing the sum of the absolute residuals around 
the plane. This method h as the advantage of being 
less affected by outliers (e.g. |j0rgensen et"aL[ll996[ ). We 

^ The SDSS-DR6 velocity dispersions are not affected by this 
bias, but they are available only for 85% of our sample. However, 
we verified that restricting the analysis to the sample with DR6 
velocity dispersions changes the FP coefficients by less than 2%. 



SDSS-UKIDSS FPs 



TABLE 1 

FP PARAMETERS IN THE r AND K BANDS FOR THE SAMPLE OF 1430 GALA 

XIES. 



RA 



DEC 



logi?e 



<A«>c 



log Re 



<A«>c,K I logo-Q 



145.34432 
147.24805 
146.81199 
146.09369 
146.46892 
146.19333 
145.68114 
145.70894 
145.48725 
145.42694 
145.44549 
145.34160 
145.19382 
146.28017 
147.30829 
146.72794 
148.85664 
147.79347 
147.74868 
148.58499 
149.12382 
149.11298 
149.18631 
148.84251 
149.11264 
149.17153 



-0.01692 
-0.03572 
-0.19005 
-0.79309 
-0.09284 
-0.03887 
-0.86722 
-0.74768 
-0.80693 
0.04954 
-0.12268 
-0.57727 
0.16887 
-0.40695 
0.15116 
-0.55688 
-0.05916 
0.12326 
0.11584 
-0.94207 
-0.39828 
-0.34883 
-0.31181 
-0.04411 
-0.47563 
-0.41298 



0.860 


21.602 


0.527 


17.046 


0.834 


21.346 


0.680 


17.437 


0.459 


19.968 


0.349 


16.319 


0.985 


21.521 


0.475 


16.216 


1.622 


23.236 


1.090 


18.009 


0.625 


20.555 


0.445 


16.439 


1.001 


21.576 


0.721 


17.350 


0.134 


18.936 


0.088 


15.653 


0.965 


21.993 


1.065 


19.167 


0.470 


19.852 


0.018 


14.473 


0.381 


19.576 


0.166 


15.351 


0.087 


18.124 


-0.152 


13.574 


0.192 


18.862 


0.137 


15.669 


0.033 


17.760 


-0.129 


13.996 


0.830 


21.046 


0.454 


16.295 


0.262 


19.223 


0.178 


15.409 


0.978 


21.643 


1.160 


19.055 


0.614 


20.337 


0.227 


15.566 


0.480 


20.362 


0.125 


15.682 


1.180 


22.575 


0.703 


17.637 


0.365 


19.682 


-0.037 


14.726 


0.685 


20.638 


0.593 


16.954 


1.291 


23.217 


0.823 


18.064 


0.469 


20.031 


0.387 


16.276 


0.572 


20.065 


0.471 


16.069 


0.359 


19.539 


0.266 


15.894 



2.282 
2.228 
2.217 
2.259 
2.317 
2.226 
2.196 
2.112 
2.180 
2.367 
2.353 
2.355 
1.930 
2.288 
2.287 
2.316 
2.226 
2.187 
2.233 
2.109 
2.324 
2.270 
2.101 
2.202 
2.398 
2.260 



adopted two fitting procedures, minimizing the residu- 
als either in log ctq or in the orthogonal direction to the 
plane. The first method is virtually independent of selec- 
tion effects in the plane of photometric parameters, while 
the orthogonal fit, adopted in most pre vious works, treats 
all th e variables symmetrically (see iLa Barbera et ahl 
[200nh . 

The FP coefficients were corrected for selection ef- 
fects through Monte-Carlo simulations. First, we gen- 
erated galaxy magnitudes according to the r-ban d lu- 
minosity function of ETGs (jBernardi et al.l[2"003af ). For 
each magnitude, we derive d log re and < m > e from the 
luminosity-size relation of IShen et al.l (|2003[ ). Values of 
log (To were assigned by using the equation of the FP, 
assuming given values of a, b, and c, and a given scatter 
in log ctq . All these quantities were chosen to match the 
observed FP. Notice that when deriving the simulated 
FP we applied the same cuts in magnitude and log ctq as 
we did for the real sample. The corrections for selection 
effects on a, b, and c, were estimated by not applying the 
cuts in magnitude and log ctq to the simulated FP. These 
corrections amount to +0.01%, +8% and +5% for the 
logCTQ fit, and -|-35% , +7% and -|-17% for th e orthogonal 
fit. As expected (e.g. lLa Barbera et al.ll200 0^. the magni- 
tude cut underestimates the coefficient a of the orthogo- 
nal fit, while for the logco fit the effect is negligible. The 
above corrections depend mainly on the scatter around 
the FP, and, because of the very similar dispersions of the 
r- and K-band FPs (see below), were applied to both the 
r- and the K-band coefficients. We notice that the above 
procedure assumes that our sample of early-type galaxies 
is magnitude complete. However, because of the match- 
ing of the initial volume-complete SDSS catalog with the 
UKIDSS database, that reduces the sample size from 
33628 to 1430 ETGs, the above assumption might not 
necessarly hold. To address this point, we retrieved ef- 



TABLE 2 

FP COEFFICIENTS. 





a 1 


b 


c 


rms 


log (To fit 

r — band 
K — band 




1.55 ±0.07 
1.51 ±0.04 


0.312 ±0.003 
0.307 ±0.003 


-9.1 ±0.1 
-8.6 ±0.1 


0.081 
0.073 


orthogonal fit 
r — band 
K — band 


1.42 ±0.05 
1.53 ±0.04 


0.305 ± 0.003 
0.308 ± 0.003 


-8.8 ±0.1 
-8.6 ±0.1 


0.064 
0.062 



fective parameters and velocity dispersions for the whole 
sample of 33628 ETGs in the SDSS catalog, and esti- 
mated how the r-band FP coefficients change between 
the whole and the UKIDSS- matched samples. We found 
the variation to be completely negligible, amounting to 
1% and to 2% for the values of a and b obtained by the 
orthogonal fit. 

The coefficients of the FP are reported in Tab. [2] In 
Fig. [U we compare the edge-on views of the r- and K- 
band planes, and in the inset we show the values of a 
and b. The a coefficients differ by only 2cr for the or- 
thogonal fit, and they almost coincide for the logcro fit, 
which is virtually unaffected by selection effects. The 
value of b ('^ 0.308) is independent of the waveband, as 
well as the scatter of the FP, which present s a tiny dif- 
ference (< 2%) for the logcro fit (see also iPahre et al.l 
|l998b) . The r-band value of a is fully consistent wit h 
that of a = 1.49 ± 0.05 found bv lBernardi et al.l (|2003bD . 
wh ile it is larger t h an th at of a = 1.24 ± 0.07 found 
bv iJorgensen et aP ()1996f ). For the K-band, the logcro 
coefficie nt is consistent with the value of a = 1.53 ± 0.08 
found bv lPahre et al.l (|1998bl), while it i s larg er than that 
ofa= 1.38±0.1 found bv lZibetti et al.l()2002f ). We notice 
that, because of the large sample size, the accuracy of our 
K-band FP coefficients is significantly higher (by 50%) 
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Fig. 1. — Short edge-on view of the r- and K-band FPs, with the photometric quantity, log Re — 6 <fi>e, being plotted as a function 
of log CO. Both quantities have been normalized by subtracting their median values, enabling a direct comparison of the r- and K-band 
projections. The median logco value amounts to 2.215dex . The FPs are represented by two lines, whose slopes are equal to the coefficients 
a of the log (TO fits. Black and grey colors correspond to the r and K bands, respectively. The average value of fe = 0.308 was adopted for 
both bands. The inset compares the r- and K-band slopes of the FP, with the concentric ellipses denoting 1 and 2cr confidence contours. 
Results of the orthogonal and logcro fits are represented by solid and dashed ellipses, respectively. 



than in previous studies. The i nvariance of the FP wit h 
waveband is in agreement with Cappella ri et al.l (120061 ). 
who found for 25 ETGs from the SAURON project the 
M/L versus L relation to have the same slope in both the 
I- and K-bands. 



5. CONSTRAINTS ON THE STELLAR POPULATIONS 

The tilt of the FP can be parameterized as a power-law 
relation between M/L and M. We assume that the stel- 
lar mass-to-light ratio of galaxies, M^/L, is a power-law: 
M^jL<x. M^ . This agrees with what was found in previ- 
ous studies for bright ETGs (see TBB04). The MjL vs. 
M r elation can then be writ ten as MjL oc M^^^ (see 
also iPrugniel fc SimieiJlTOQ^ . The quantity /3 is related 
to the ratio of stellar to total mass as M^/M oc M~^, 
describing the variation of wavelength independent prop- 
erties with mass (e.g. dark matter content). The quan- 
tity /3* depends on how stellar population parameters 
vary with mass. Considering only the effects of age and 



metallicity, for each waveband, we can write: 
Si\ogAU/L) 13* 5l\ogt) , 5{\ogZ) 



5{\QgM,) 



1-/3 



Ct 



<5(logM*) 



cz 



5{\ogNU 



\) 



where the quantities 5{\ogt) and 6{\ogZ) are the log- 
arithmic differences of age and metallicity defined be- 
tween more and less massive galaxies (per decade in stel- 

lar mass), while ct = gjog/ and cz = qio^z ^"^^ 
the partial derivatives of logM^/L (in a given band) with 
respect to t and Z. Writing Eq. [l] for r and K bands, 
we obtain two independent constraints on (5(logi) and 
S{\ogZ). Then, expressing the values of P* in the r and 
K bands (/3* and (3^) as a function of the correspond- 
ing coefficients of the FP, we can estimate 6{\ogt) and 
6{\ogZ). We introduce the parameter / = I3'^I{I3 + (3^) 
which defines the fraction of the K-band slope of the 
M/L vs. M relation (i.e. the fraction of the K-band 
tilt) due to stellar population effects. We note that / can 
vary between and 1. For / = 0, the K-band tilt is 
independent of stellar populations (/3^ = 0), while for 
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Fig. 2. — Differences of age, 5(logt), and metallicty, 5(logZ), per decade in mass between more and less massive galaxies, as derived 
from the FP slopes in the r and K bands, by adopting the mean values of a of the orthogonal and logcro fits (see the text). Red crosses and 
blue circles correspond to the cases of homology assumption with / = and / = 1, respectively, while the green squares correspond to the 
case where non-homology is taken into account (see Sec. 5). For each symbol type, different points mark the values of 5{logt) and <5(log Z) 
obtained by shifting a^ and ax according to the corresponding uncertainties. For each case, the median values, St and <5Z, of S{logt) and 
(5(log Z) are plotted as black symbols and are reported in the lower-left part of the plot. The effect of changing / on the median values of 
<5(log t) and (5(log Z) is shown by the black solid and dashed lines, for the BrC03 and the Vaz96 models, respectively. The horizontal dotted 
lines mark the values of <5Z obtained from previous studies of the color-magnitude relation (KOD98) and absorption line indices (THO05). 



/ = 1 the tilt is entirely explained by stellar population 
effects {(3 = 0). Under the assumption of homology, the 
slope of the M/L vs. M relation can be written as 
{2 — a)/ {2 + a). With the above notation, the following 
relations hold: 

/? + /3; = (1 - /)//-/?^ + A* = (2 - a,)/{2 + a,), (2a) 
P + 0K= Pk/f = (2 - flK)/(2 + ak). {2b) 

For different /'s and using the values of a in Tab. [21 we 
computed (3* and P^ from Eqs. [H and then, inverting 
Eq. [T]for both bands, we derived 6{\ogt) and i5(log^). 
The values of c^ and c-z were estimated using simple stel- 
lar population models, with s olar metallicity and an age 
of 12 Gyr ^, using both the iBruzual fc CharlotI (|2003|) 
(from now BrC03) and the updated IVazdekig ()1996D 



^ Computing Ct and cz by varying t between 5 and 12Gyr, and 
Z from 0.5Zq to 2.0Zq, we found that the change of the median 
values of 5(logt) and 5(logZ) are negligible, amounting to less 
than O.Oldex and O.OOSdex, respectively. 



(from now Vaz96) models. We adopted a Scalo IMF and 
a Salpeter IMF for the BrC03 and Vaz96 models, respec- 
tively. Fig.[2]shows (5(logZ) versus (5(logi) obtained for 
/ = and / = 1, as well as the mean values of S{\ogt) 
and S{\ogZ) as a function of /. The scatter seen in the 
plot reflects the uncertainties listed in Tab. [5] for the a^- 
and Ok coefficients. The figure shows that if the NIR 
tilt of the FP is not caused by stellar population effects 
(/ = 0), more massive galaxies have to be more metal 
rich than less massive ones {S{\ogZ)> 0), with galaxy 
ages being remarkably homogeneous {St/t '^1%). As / 
increases, we see that S{logZ) decreases, while S{\ogt) 
becomes larger. Specifically, for / '^ 1, more massive 
galaxies are much older {St/t ^ 50%) and less metal rich 
than low mass systems. 

6. DISCUSSION 

This work presents the waveband dependence of the 
FP by comparing the optical and NIR FPs for a large 
sample of galaxies, with homogeneous measurements of 
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structural parameters and velocity dispersions. This is 
allowed, for the first time, thanks to the availability of 
both SDSS r-band photometry and spectroscopy, and 
UKIDSS K-band photometry for the same sample of 
ETGs. Such a dataset, together with the use of the same 
fitting procedure in both bands, makes our study virtu- 
ally free from any methodological effect on the waveband 
dependence of the FP. Our analysis shows that the FP 
does not change significantly from the optical to the NIR, 
bringing interesting questions about the nature of the se- 
quence of ETGs. 

In Sec. 5, we have shown how the r- and K-band FPs 
constrain the variation of the stellar population prop- 
erties (age and metallicity) with stellar mass, and how 
such a constraint is strongly dependent on the fraction, 
/, of the FP tilt resulting from stellar population ef- 
fects. Previous studies of the color-magnitude (CM) re- 
lation and of line-strength indices of ETGs might help 
us to s olve this dependency, by deriving the proper value 
for /. iKodama et al.l ()1998f ) showed that the little red- 
shift evolution of the CM relation implies (i) all the 
(luminous) ETGs to be equally old and (ii) more mas- 
sive galaxies to be more metal rich, with a metallic- 
ity change of S{\ogZ) ^ 0.22dex per decade in stellar 
m ass. This fin ding is qualitatively consistent with that 
of iThomas et a l. (2005), who found absorption- line in- 
dices consistent with a metallicity change of 5{\ogZ) ~ 
0.12dex per mass decade. They also found evidence 
for an age gradient along the sequence of ETGs, with 
^(logi) = -h0.05 ± 0.07 (see their eq.3). Fig. [2] com- 
pares these values of 6{logZ) with those derived by the 
optical-NIR FP. For the BrC03 (Vaz96) model, the max- 
imum value of 6(\og Z), which is consistent with the FP, 
amounts to 0.06(0.09) ±0.04(0.04)dex for / = 0. This 
value is 4cr (3.7(t) lower than that derived by the CM re- 
lat ion, while it is only 1 .5a (O.Scr) lower than that found 
bv lThomas et al.l (|2005f l. As / increases, the FP requires 
the value of S{log Z) to decrease, making the above dif- 
ferences even larger. Therefore, reconciling previous es- 
timates of 6{logZ) with our results leads to a scenario 
where / = 0, which means that the FP tilt is not driven 
by stellar populations. We have to remark, however, 
that this interpretation is troublesome, since galaxy col- 
ors and line indices are always measured within a given 
fixed aperture, and the presence of internal population 
gradients in galaxies can significant ly affect the inf erred 
values of (5(logZ) and 6{logt) (e.g. iScode ggio 20011). 

Further constraints come from previous works ad- 
dressing the origin of the FP tilt itself. Performing 
a detailed dyna mical analysis of twenty- five galaxies, 
ICappellari et al.l ([2006) derived M/L ratios consistent 
with those obtained from the virial theorem in the as- 
sumption of homology, concluding that structural and 
orbital non-hom ology have a negligib le role in the tilt of 
the FP (see also IZaritskv et ai]l2008[ ). In support to this 
view, they also showed that the variation of the dynam- 
ical M/L is correlated with the H^ line-strength, thus 
ascribing most of the tilt to stellar population (age) ef- 
fects. However, as the authors notice, this result strictly 
applies to their measurement of the velocity dispersion 
as the average over an aperture of radius equal to r^, 
a fact that alone might compensate part of the dynam- 
ical non-homology. Moreover, most of the galaxies in 
their sample (68%) are fast rotators, while five of them 



(20%) have low velocity dispersio n (a =:60-85 km s ^). 
As fou nd by IZaritskv et aH i|2006D and iD'Onofrio eFaD 
()2008[) . bright and faint spheroids have different FPs, 
with the tilt becoming larger for galaxies having low 
velocity dispersion. Hen ce, the different selection of 
our sample and that of ICappellari et al.l (l2006l) pre - 
vents a straightforward comparison. iBolton et al.l ()2007f ) 
showed that by replacing mean surface brightness with 
mass density, the FP relation closely approaches the 
virial theorem expectation, implying that most of the 
tilt is caused by a variation of dark matter content 
with galaxy mass. However, the uncertainties on their 
FP coefficients and the possible biases introduced by 
the gravitational-lens selection do not definitively ex- 
clude the contribution of non-ho mology t o the FP 
tilt. Our result agrees with Bolton et al.l ()2007| ) re- 
garding the minor role played b y stellar populat i ons on 
the tilt. TBB04, agreeing w ith iBusarello et~all ()1997f ) 
and lGraham fc CoUessI ()1997| ). found that structural and 
dynamical non-homology can account for more than two- 
thirds of the FP tilt, with the remaining part being ex- 
plained by stellar population effects. In particular, re- 
stricting the analysis to the magnitude-complete subsam- 
ple, they found that the contribution of stellar popula- 
tions to the tilt becomes negligible. Notice also that their 
FPs are derived from different sources, with significantly 
different coefficients in the optical and NIR wavebands. 

To understand how our results may be affected by the 
assumption of homology, we followed an approach similar 
to that of TBB04, using spherical, isotropic, non rotat- 
ing, one- component models of ET Gs following the Sersic 
law (see iLa Barbera et al.l l2005t ). For each galaxy, we 
considered the model with the corresponding Sersic in- 
dex in the r-band, and used that model to correct the 
central velocity dispersion to the quantity at (defined as 
the square-root of the total specific kinetic energy) , the 
effective radius to the gravitational radius, Vg, and to 
calculate the mean surface brightness within the grav- 
itational radius, < /i >g. Applying the orthogonal fit, 
we obtain the following equation of the FP in the K 
band: logrg oc (2.3 ± 0.2) logcr* + (0.4 ± 0.02) < ^ >g, 
which is remarkably close to the virial theorem expecta- 
tion, implying that non-homology may account for the 
entire tilt. To explore how this result would affect stellar 
population properties, we normalized the r- and K-band 
FP coefficients in such a way to match the virial theo- 
rem expectation (a = 2.0 and b — 0.4) in the K band. 
Then, we derived the corresponding values of S{logZ) 
and d{logt) (see Sec. [5]). As we can see in Fig. [2l ac- 
counting for non-homology leads to the same S{\ogZ) 
and 6(logt) values as those derived under the assump- 
tion of homology if the tilt of the FP is not due to stellar 
populations (i.e. / = 0). As discussed above, this is 
consistent with what expected from the color-magnitude 
relation and absorption-line indices of ETGs. 

In summary, our analysis suggests a consistent picture 
where (i) the tilt of the FP does not originate from stellar 
population effects but is due to other effects, such as non- 
homology; (u) the SDSS-UKIDSS FPs require more mas- 
sive galaxies to be mildly more metal rich than less mas- 
sive systems, and to have extremely synchronized ages, 
with the age variation per mass decade being smaller 
than few percent. 
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